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Application of modified simple 
equation method to Burgers, Huxley 
and Burgers-Huxley equations 


Z. Ayati*, M. Moradi and M. Mirzazadeh 


Abstract 


In this paper, modified simple equation method has been applied to ob- 
tain generalized solutions of Burgers, Huxley equations and combined forms 
of these equations. The new exact solutions of these equations have been 
obtained. It has been shown that the proposed method provides a very effec- 
tive, and powerful mathematical tool for solving nonlinear partial differential 
equations. 


Keywords: Modified simple equation method; Burgers equation; Huxley 
equation; Burger-Huxley equation. 


1 Introduction 


Mathematical modeling of many real phenomena leads to a non-linear par- 
tial differential equations in various fields of sciences and engineering. Many 
powerful methods have been presented for solving PDEs so far, such as 
tanh-function method [19] and [28], sine-cosine method [29], Homotopy 
Analysis method [17], Homotopy perturbation method [6], variational itera- 
tion method [9] and [10], Adomian decomposition method [1], Exp-function 
method [1], [11], [36] and [37], simplest equation method [7] and [4], and 
many others. Most recently, a modification of simplest equation method 
(MSE method) has been developed to obtain solutions of various nonlinear 


*Corresponding author 

Received 3 June 2014; revised 7 April 2015; accepted 17 June 2015 

Z. Ayati 

Department of Engineering Sciences, Faculty of Technology and Engineering East of 
Guilan, University of Guilan, Rudsar-Vajargah, Iran. e-mail: zainab.ayati@guilan.ac.ir 
; Ayati.zainabQ@gmail.com 


M. Moradi 
Department of Engineering Sciences, Faculty of Technology and Engineering East of 
Guilan, University of Guilan, Rudsar-Vajargah, Iran. email: 


M. Mirzazadeh 
Department of Engineering Sciences, Faculty of Technology and Engineering East of 
Guilan, University of Guilan, Rudsar-Vajargah, Iran. email: 


59 


60 Z. Ayati, M. Moradi and M. Mirzazadeh 


evolution equations [14], [15], [21], [31], [82], [33] and [34]. The present paper 
is motivated by the desire to extend the MSE method to obtain generalized 
solutions of Burgers, Huxley, and Burgers-Huxley. The procedure of this 
method, by the help of Matlab, Maple or any mathematical package, is of 
utter simplicity. 


2 The MSE method 


Consider a nonlinear partial equation in two independent variables, say x 
and t , in the form of 


P(U, Ut, Un, Ute, Ure, aa) = 0, (1) 


where u = u(x,t) an unknown function, P is a polynomial in u = u(z, t) 
and its various partial derivatives, in which the highest order derivatives and 
nonlinear terms are involved. This method consists of the following steps. 


Step 1. Using the transformation 
€=2+ut, (2) 


where w is constant, we can rewrite equation (1) as a following nonlinear 
ODE: 


Q(u,u',u”’,...) =0. (3) 


Where the superscripts denote the derivatives with respect to € . 


Step 2. Suppose that the solution of equation (2) can be expressed as 


follows “ 
ul) = al Fey. (4) 
1=0 


Where a; are constants to be determined later, with am 4 0 and F(€) is an 
unknown function to be determined later. 


Step 3. The positive integer m can be determined by considering the homo- 
geneous balance of the highest order derivatives and highest order nonlinear 
appearing in equation (2). 


Step 4. Calculating all necessary derivatives u’,u”’,u'”’,... , and substituting 
F'() 
; F(§) 
Equating the coefficients of same power of F~*(£) to zero gives a system of 


equations which can be used to solve for determining unknown constants,F (€) 


equation (3) into equation (2) yields a polynomial of and its derivatives. 
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and F’(€). By substituting obtained results into equation (3), solutions of 
the equation (1) can be obtained. 


3 Application of the MSE method 


In this section, the modified simple equation method has been applied to 
obtain generalized solutions of Burgers, Huxley, and Burgers-Huxley. 


3.1 Application MSE method to Burgers equation 


The Burgers equation is a nonlinear partial differential equation of second 
order of the form 
Ut + UWUy = VUge- (5) 


Where v is the viscosity coefficient [2], [22], [23] and [27]. Many problems 
can be modeled by the Burgers‘ equation. This equation is one of the very 
few nonlinear partial differential equations which can be solved exactly for 
the restricted set of initial function. The study of the general properties of 
the equation has drawn considerable attention due to its place of application 
in some fields such as gas dynamics, heat conduction, elastically, etc. 
To apply MSE method on equation (4), lets introduce a variable € , defined 
as 

€=x2- ut. (6) 


So, equation (4) turns to the following system of ordinary different equation, 

—wu +uu' = vu". (7) 
Where w is constant to be determined. Integrating (7) and considering the 
integral constant to be zero, we obtain 


—wut ze =vu'. (8) 


Suppose that the solution of ODE equation (8) can be expressed by a poly- 
nomial in E as shown in (3). Balancing the terms u? and u’ in equation (8), 
yields to m = 1. So we can write (3) as the following simple form 


u(6) = an Fee} + aos # 0. (9) 
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So 
wu! =ai(— — (=)’). (10) 


Substituting (9) and (10) into equation (8) and equating each coefficient of 
F-*(€) , (« =0,1,2) to zero, we derive 


1 
—wag + 5% = 0, (11) 
(—w + a9) F’ —vF" =0, (12) 
1 
(541 + va1)(F’)? = 0. (13) 


By solving equations (11) and (13), the following results will be obtained 


do = 0, 2w,a, = —2v. 


Case 1. when equation (12) turns to 
wk’ +vF"=0. 


So 7 
F’ = Aew§, (14) 
Where A is a arbitrary constant. Integrating (13) with respect €, F(€) will 
be obtained as fallows 
—Av 
F — 
w 


ev § + B, 


where B is a constant of integration. Now, the exact solution of equation 
(4) has the form 
—QvAer @-¥) 
Av 04 (x—ut) a B S 


Ww 


ui(a,t) = 


Case 2. when ap = 2w, equation (12) yields to 


wk’ —vF" =0. 
So 
F' = Aes, 
and A 
PEE eee, 
w 


Now, the exact solution of equation (4) has the form 
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Qwy Aer (wt) - 2w?B 
Aver @-¥) 4 wB Aver (e-vt) + wB’ 


ug(a,t) = 2w 


Note that all obtained solutions have been checked with maple 13 by 
putting into the original equation and found correct. 


3.2 Application MSE method to Huxley equation 


Now we will bring to bear the MSE method to obtain exact solution to the 
Huxley equation [3], [7], [8], [12], [13], [18], [25] and [26] in the following form 


Ut = Ure + u(k — u)(u— 1). (15) 


The Huxley equation is an evolution equation that describes the nerve prop- 
agation in biology from which molecular CB properties can be calculated. 
It also gives a phenomenological description of the behaviour of the myosin 
heads II. This equation has many fascinating phenomena such as bursting 
oscillation [3], interspike [18], bifurcation, and chaos [35]. A generalized ex- 
act solution can gain an insight into these phenomena. There is no universal 
method for nonlinear equations. In this part, the exact solution will be ob- 
tained by the MSE method. 
By considering (6), equation (15) turns to the following ordinary differ- 
ential equation, 
wu +u"+u(k—u)(u-—1) =0. (16) 


Balancing the terms u” and u® in equation (16), yields to m = 1. So we 
can rewrite (3) as the following simple form 


FQ) 
‘FE 


Now by substituting 19) into equation (16) and equating each coefficient of 
F~‘(€) ,(i = 0,1,2,3) to zero, the following result will be obtained 


u(€) =a + ao, a1 #0. (17) 


ag — (k +1)a2 + kag = 0, (18) 

a, F’" + wa, F" + (2(k + 1)agai — 3aga1 — kay) F’ = 0, (19) 
—3a,F'F" + ((k + 1)a, — wa, — 3aqa?)(F’)? = 0, (20) 
(2a, — a3)(F’)? = 0. (21) 


Solving equations (15) and (22), we drive 
ag = 0, 1, k, 
a= +2. 
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Case 1. when ao = 0, equations (16) and (17) yield 
F" + wF" —kF' =0, (22) 


3F" + (w—k-1)F' =0. (23) 
By substituting equation (23) into (22), we obtain 


3k 
pes eS, 
ee w—-k— iD 
So 
Fl = Ae. (24) 


Where a = w+ a and A is a arbitrary constant. Therefore, we have 


A 
F'= sae i +B. (25) 


where A and B are arbitrary constants. By substituting (25) into equations 
(22) and (23), we get 


k+1. 3 
= + k? k+1,B=0. 
7 7 6 ; 0 


Thus, (25) can be rewritten as follows 


P= A at. (26) 
a 


Integrating (26) with respect € , F(&) will be obtained as follows 


A 
F = sens + C,; 
a 


where C is a constant of integration. Substituting the value of F and F” into 
equation (19), the following exact solution of equation (16) has been obtained 


+,/2 Aae~% 
Ae~% + a2C" 


ui(a,t) = 


Case 2. when ag = 1, equations (16) and (17) turns to 
Fl" +wk" +(k-1)F"=0, (27) 


3F" + ((w—k—-143V2)F’ =0. (28) 
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By substituting equation (8) into (27), we obtain 


3(k — 1) 
w—k—-143V2 


FM + (w )F" =0. 


So 
Be Ae- (29) 


ee and A is a arbitrary constant. Therefore, we have 


where a@ = w— ek eae 


A 
Fi =——e-% +B, (30) 
a 


where A and B are arbitrary constants. By substituting (30) into equations 
(27) and (28), we get 


Rea 9 3 
—_ 2 _ 
w= ava sh 6k — 6kV/2 + 27 — 6V2, B =0, 


or 


k+1 
w= SVD ye 6k + 6kV2 + 27 + 62, B =0. 


Thus, (30) can be rewritten as follows 


Fo =——e%, (31) 


A 
a 
Integrating (31) with respect , will be obtained as follows 


A 
F = ses + C, 
a 


where C is a constant of integration. Substituting the value of F and F” into 
equation (19), the following exact solution of equation (16) has been obtained 


+/2 Aae~% 


Case 3. when ao = k, equations (16) and (17) turn to 


F" 4+-wF"” +k(1—k)F’ =0, 


3F" + ((w—14+k(-1+3V2))F’ =0. 
By using the same method applied in case 1, the following solution will 


be obtained 
4 +/2Aae~% 
Ae~% + a2C’ 


u3(a,t) =k 
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3k(1—k) 
w—1+k(-1+3V2)’ 


w= ah 3v/2) a oy 2k(-1+ 3V2) +k(-1+ 3/2)? — 8k? + 8k. 


3.3 Application MSE method to Burgers-Huxley 
equation 


The analysis presented in this part is based on the generalized nonlinear 
Burgers-Huxley equation, 


Ut = Ure + UWUe + u(k — u)(u— 1), (32) 


which models the interaction between reaction mechanisms, convection ef- 
fects and diffusion transports [20], and some special cases of the equation, 
which usually appear in mathematical modelling of some real world phenom- 
ena. It also gives a phenomenological description of the behaviour of the 
myosin heads IT [30] and Fitzhugh-Nagoma equation, an important nonlin- 
ear reaction-diffusion equation used in circuit theory, biology and population 
genetics [5]. 

By considering (6), equation (32) turns to the following ordinary differ- 
ential equation, 


wu tu" +uu' +u(k —u)(u—1) =0. (33) 


Balancing the terms wu” and u? in Eq. (33), yields to m = 1. So we can 
rewrite (3) as the following simple form 
Uy 


u(€) = a5) + ao, #0. (34) 


Now by substituting (34) into equation (33) and equating each coefficient of 
F~‘(€), (i =0,1,2,3) to zero, the following result will be obtained 


aa — (k + 1)a? + kag = 0, (35) 


a, F"” + (way + ajay) F” + (2(k + 1)agay — 3a3a, —ka,)F’=0, — (36) 
(—3a, + a3) FF” + ((k + 1)a, — way, — aga, — 3aga7)(F’)? = 0, (37) 
(2a; — a? — a3)(F’)? = 0. (38) 


Solving equation (35) and (38), we get 
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ag = 0, 1, k, 


ay = 1, —2. 


Case 1. when ag = 0 and a, = 1 equations (36) and (37) yield 
F" 4 wF" —kF’ =0, (39) 


2F" +(w—k—1)F’ =0. (40) 


By substituting equation (40) into (39), we obtain 


2k 
pl’ | fF” = 0. 
w w—k—- i) 
So 
F” = Ae~% (41) 
where a= w+ — 2k) and A is a arbitrary constant. Therefore, we have 


A 
F'= “eee +B, (42) 


where A and B are arbitrary constants. By substituting (42) into equations 
(39) and (40), we get 
w=4+(k—1),B=0, 


So 
a=-—1l,-k. 


Thus, (42) can be rewritten as follows 


A Ok 


Fo = ——e 43 
~e (43) 


Integrating (43) with respect € , F(&) will be obtained as follows 


A 
F' = sens +C, 
Qa 


where C is a constant of integration. Substituting the value of F and F’ 
into equation (34), the following exact solutions of equation (33) has been 
obtained 


Aet—(k-Dt 
u(x,t) = Act—(k-Dt 4 ©’ 
AK ck(e+(k-1)t) 
ua(@, t) = 


= Aek(a+(k-Lt) 4 6200 
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Case 2. when ap = 0, and a, = —2 equations (36) and (37) yield 


F" +wF" —kF' =0, (44) 


5F" + (w—k-—1)F"” =0. (45) 
By substituting equation (45) into (44), we obtain 


5k) 
w—k—-1 


Fl + (w4 le 0s 


By using the same method applied in case 1, the following solution will be 
obtained 


2Aae~% 
walt) = 7 omok + GBC 
where 
5k 
a=wt Rees ae 
and 3 5 
w= 3(k t 1)+ vk? 14k +1. 


Case 3. when ag = 1 and a, = 1, equations (31) and (32) turn to 


Fl +(w+1)F" 4+ (k-1)F’ =0, (46) 


2F" + (w+3-—k)F’ =0. (47) 


By substituting equation (47) into (46), we obtain 


2(k —1) 
Fl’ I | 1 FF" = 0 
( w+3— Be 
So 
F" = Ae~%, 
where a = w+1-— 2h) and A is an arbitrary constant. Therefore, we have 


A 
F'= a +B, (48) 


where A and B are arbitrary constants. 
By substituting (48) into equations (46) and (47), we get 


w=-1tk,B=0. 
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Thus, the following exact solutions of equation (33) has been obtained. 


Ae~2t+U—kyt 
us (x, t) =1 Ae-#t+G-h)t 4 ©? 


A(1 oa k)e~A-k) e+ A-+k)#) 
U6 (2, t) =1 Ae (1-k)(a@+(+k)t) 4 (1 - k2C" 


Case 4. when ag = 1 and a; = —2, by using the same method applied 
in case 1, the following solution will be obtained. 


2Aae~% 


Sy Se 
urs(z,t) = Ae~% + a2C 
Where 5(k — 1) 
Be Ey ara ae 
a 13 + 3k 5 
Pea LT tes pony SRT PTO 
8 8 
Case 5. when ao = k and a; = 1, equations (31) and (32) turn to 
FY + (wt+k)F" +k(1—k)F’ =0, (49) 
2F" + (w+3k—1)F" =0. (50) 


By substituting equation (50) into (49), we obtain 


2k(1 —k) 
w+3k—-1 


FU +(wtk VE =O. 


So 
F" = Ae~%, 


where a = w+k-— 2k) and A is an arbitrary constant. Therefore, we 


have 


A 
Fi => ——e-“% + B, (51) 
(a7 


where A and B are arbitrary constants. By substituting (51) into equations 


(49) and (50), we get 
w=—-k+1,B=0. 


Thus, the following exact solutions of equation (32) has been obtained 


Ake—k(@+(k- 10) 
ug(a,t) = k Ac k@HE-Di) 4 k20" 
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A(k — Pere Yeteeyy) 
Aen (k-1)(a@+(k+1)t) + (k —s 1)2C’ 


uio(a, t) = k 


Case 6. when ap = k and a, = —2, by using the same method applied 
in case 5, the following solution will be obtained. 


2Aae~ 
ui1,12(#,t) =k+ Aen 4p 08O" 
Where 5k(1 — k) 
sea Fa TS 
and 
w= ate t1)+ 2 33 + 30k +1. 


4 Conclusion 


In this paper, modified simple equation method has been applied to obtain 
the generalized solutions of some nonlinear partial differential equation. The 
results show that modified simple equation method is a powerful tool for 
obtaining the exact solutions of nonlinear differential equations. It may be 
concluded that, the method can be easily extended to all kinds of nonlinear 
equations. The advantage of this method over other methods is that in most 
methods applied for the exact solution of partial differential equations such 
as Exp-function method, e - expansion method, tanh-function method, and 
so on, the solution is presented in terms of some pre-defined functions, but in 
the MSE method, F(€) is not pre-defined or not a solution of any pre-defined 
equation. Therefore, some new solutions might be found by this method.The 
computations associated in this work were performed by Maple 13. 
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